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Abstract 

We provide a brief survey of quantum statistical characterisations of order, disorder 
and coherence in systems of many degree of freedom. Here, order and coherence are 
described in terms of symmetry breakdown, while disorder is described in terms of entropy 
CN) ■ and algorithmic complexity, whose interconnection has been recently extended from the 
classical to the quantum domain. We see that, in the present physical context, the concepts 
'. of order and disorder are not mutually antithetical but bear an interrelationship similar to 
^p H ' that between signals and noise. 

•*-> . 

Sh 1 1. Introduction. 

This article is designed to provide a succinct account of the prevailing quantum sta- 
tistical pictures of order and disorder in systems with many degrees of freedom. In this 
context, order essentially signifies organisation of the microscopic components of such sys- 
tems to produce macroscopic fields, or signals, as exemplified by the polarisation of a 
^ ■ ferromagnet; while disorder amounts to randomness. 

Here we formulate mathematical pictures of order and disorder within the framework 
of operator algebraic statistical mechanics [Em, Th, Sel] , which provides a natural setting 
for their descriptions. We start in Sec. 2 with a brief sketch of the structure of algebraic 
quantum theory. We then pass on, in Sec. 3, to both the probabilistic formulation of 
^ disorder in terms of Von Neumann's entropy [VN] and the intrinsic description thereof by 

Kolmogorov's algorithmic complexity [Ko]. In particular, we discuss Brudno's theorem 
[Br] and its recent quantum generalisation [BKMSS], which shows that these two char- 
acterisations of disorder essentially yield the same picture. In Sec. 4, we formulate the 
concept of order due to symmetry breakdown. In Sec. 5, we refine this formulation of 
order to an extreme version thereof, namely coherence, in the sense proposed by Glauber 
[Gl]. We provide some concrete examples both of order, in Sec. 4, and of coherence in 
Sec. 5. We conclude in Sec. 6 with some further brief observations about order, disorder 
and coherence, and discuss the need to widen the concept of order to the description of 
organisational structures that are not covered by existing theories. 

2. The Operator Algebraic Framework 

We employ the standard operator algebraic description [Em, Th, Sel] of a quantum 



O 



i 



mechanical system, E, as a triple (A,S,a) representing its observables, states and dy- 
namics, respectively. Specifically, A is a C*-algebra, whose self-adjoint elements represent 
the bounded observables of E, and a is a homomorphism of either the additive group R 
into the automorphisms of A or of the semigroup R + into completely positive, identity 
preserving, linear contractions of this algebra, according to whether E is a conservative 
system or an open Markovian dissipative one. The state space S is a norm closed, convex 
subset of the positive, normalised, linear functionals on A that is stable under the action 
of the dual of a. We shall denote the expectation value of A (&A) for the state p by p(A) 
or, equivalently, (p; A). The pure states are the extremal elements of iS.Thus, the model is 
specified by the structures of A, S and a. We note that this generic model also covers the 
case of classical mechanical systems, which are distinguished by the condition that A is 
abelian. In this case, by the Gelfand isomorphism, A is the algebra of continuous functions 
on a compact space K, S is a set of probability measures on K and the transformations 
at are implemented by transformations Tt of K. Here K is the 'phase space' of the model. 

The Finite System Model [VN]. For this, A is the VF*-algebra of bounded operators 
in a Hilbert space TL and S is the set of normal, i.e. ultraweakly continuous, states, p, on 
A: these correspond to density matrices, denoted by the same symbol, according to the 
formula p(A)=Tr{pA). Thus S is a convex set and its extremal elements, representing the 
pure states, are those whose density matrices are one-dimensional projectors. In the case 
where E is conservative, its dynamical automorphisms, at, are unitary transformations 
A^ex.p(iHt)Aex.p(— iHt) of A, where H is the Hamiltonian operator of E in units for 
which h = 1. In the case where the system is dissipative and its dynamical semigroup a is 
strongly continuous, its generator L takes the following form [Li] . 

LA = i[H,A}_ + J2 r (VrAV r - ±[V?V r ,A]+), (2.1) 

where H (= if*), V r and ^2 r V*V r belong to A and [., .] T denote commutator and anti- 
commutator, respectively. 

The Infinite System Model [Em, Th, Sel]. This represents a system, E, of particles 
that occupies an infinitely extended space X, which we take to be either a Euclidean 
continuum, R d , or a lattice, Z d . We denote by C the set of all bounded open regions of 
X, and, for each A in £, we construct a VF*-algebra, Aa, of observables that is just that 
of a system, Ea, of particles of the given species confined to A. These local algebras are 
constructed so as to satisfy the natural demands that Aa is isotonic with respect to A 
and that Aa and Aa' intercommute if A and A' are disjoint. It follows from the isotony 
property that Ac '■= Ua6£^A' ^ s well-defined normed *-algebra. We designate its norm 
completion, A, to be the C*-algebra of the bounded observables of E. We assume that 
this algebra is equipped with a representation, 7, of the space translation group X in its 
automorphisms, which satsifies the covariance condition that , j(x)Aa=Aa+x- 

We assume that the state space, S, is a convex set of positive, normalised, linear func- 
tionals on A, whose restrictions to the local algebras Aa are normal: the local normality 
condition serves to exclude the possibility of finding an infinity of particles in a bounded 
spatial region [DDR]. 
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The dynamics of E is formulated as a natural infinite volume limit of that of the finite 
system Ea- In the conservative case, this latter dynamics is governed by the form of the 
Hamiltonian operator, Ha, affiliated* to Ax- Thus, if A is an element of Ac and therefore 
of A\ for A sufficiently large, exp(iHAt)Aexp(— iH^t) is its evolute at time t with respect 
to the dynamics of Sa- In the simplest cases, such as that of lattice systems with short 
range interactions [St, Ro], this coverges in norm to a definite limit as A increases to X 
over a sequence of suitably regular regions and thus yields a definition of the dynamical 
automorphisms a by the formula 

a t A = norm : lim. A ^x^p{^H\t)Aexp(-iH A t) V AeA c , teR. (2.2) 

More generally, when the convergence condition for this formula is not fulfilled, S has to be 
formulated so as to comprise just those states that support a limit dynamics represented 
by a weaker form of Eq. (2.2) [Sel, 2]. In the case where E is an open dissipative system, 
its dynamical semigroup is similarly defined as a limit of that of the corresponding finite 
system Ea- 

Thus the model of E, is represented by the quadruple (^4, 5, 7, a). The subset Sx of 
S comprising the space translationally invariant states is manifestly convex and we denote 
by £(Sx) the set of its extremal elements. These are termed the spatially ergodic states. 

Affiliated Quantum Fields. Identifying the algebra A with any faithful represen- 
tation thereof, a quantum field, £(x), of the model is defined to be a distribution valued 
operator that is covariant with respect to space translations and whose integral against a 
test function f(x) with compact support is affiliated to the local algebras A A for which 
ADsupp(/). 

Explicit Constructions. The local algebras Aa, the space translation group 7 and 
the local Hamiltonians Ha, on which the model of E is based, are constructed as follows 
(cf. [St, Ro, HHW] or the general treatments [Em, Th, Sel]). 

In the case where E is a system of particles, e.g. Pauli spins, on a lattice X = Z d , 
we assume that the algebra of observables, Aq, of each particle is that of the operators 
in a finite dimensional Hilbert space, Hq. We take the local algebra Aa, for Ae£, to be 
the tensor product ®x&aA x of copies A x of Aq attached to the respective sites x in A; 
and, for AcA', we identify A (<eAa) with A®I a >\a (eAa 1 )- Under this identification, the 
algebras {^4a|Ag£} satisfy the conditions of isotony and local commutativity and thus 
permit the above definitions of Ac and A. Further, denoting by a x the copy in A x of 
the element ao of Aq, we define the space translational automorphism group 7 by the 
formula r ){x)a x i = a x+x >. Thus, the local algebras {Aa} transform covariantly under this 
group. The local Hamiltonian Ha is the element of Aa representing the interaction energy 
involving only the particles in A. 

In the case where E is a system of particles of one species in the Euclidean continuum 
X = R d , we formulate its observables in second quantisation, as expressed in terms of a 

* A possibly unbounded operator Q in the representation space of a IU*-algebra B is 
said to be affiliated to B if it commutes with B' ', the commutant of B. 
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quantised scalar or spinor field ip, according to whether the particles are bosons or fermions. 
In either case, ip is a distribution valued operator in Fock space Tip. Its algebraic properties 
are governed by the canonical commutation or anticommutation relations according to 
whether E is composed of bosons or fermions. We define the local Hilbert space Ha to be 
the subspace ofH,F generated by application to the Fock vacuum of the polynomials in the 
smeared fields obtained by integrating tp* against P(A)-class test functions and we define 
the local algebra Aa to be that of the bounded operators in H,a- The space translational 
automorphism group 7 is defined by the canonical formula , y(x) , ifj(x')= , ifj(x + x') . The 
local Hamiltonian Ha is just that of Ea and we assume that it, and consequently also the 
automorphisms at, is invariant under the gauge automorphisms ip(x)—>ip(x)exp(ic), with 
c real and constant. 

Equilibrium States of Conservative Systems. The equilibrium states of a conser- 
vative system E at inverse temperature (3 are characterised by the Kubo-Martin-Schwinger 
(KMS) condition, namely (cf. [HHW, Em, Th, Sel]) 

(p; [a t A]B) = (p; Ba t+i/3 A) V A, BeA, teK. (2.3) 

This represents various conditions of dynamical and thermodynamical stability [Sel] and, 
in the case of an infinite system, it automatically ensures that p is locally normal [TW]. 
In general, it follows from Eq. (2.3) that its equilibrium (KMS) states at the inverse 
temperature (3 comprise a convex set, which we denote by Sp. In the case where E is a finite 
system, Sp consists of just the canonical state with density matrix exp(— / 9if)/Tr(Idem). 
By contrast, for an infinite system, Sp is a Choquet simplex, which may contain more 
than one element and whose decomposition into extremals is just the central one [Ru, 
EKV]. Thus the set £{Sp), of its extremals consists of primary states and may naturally 
be interpreted as comprising the pure equilibrium phases [EKV] of the system. Moreover, 
as they are primary, they enjoy the clustering property that [Ru] 

lim| xHoo [(p; A>y(x)B) - (p; A)(p; 7 (x)S)] = V A, BeA. (2.4) 

Open Dissipative Systems. The situation is different for these systems since they 
carry no natural counterpart of the KMS states. In particular, the model (A, S, 7, a) does 
not necessarily have any stationary, dynamically stable primary states, which could be the 
counterparts of the pure phase equilibrium states of conservative systems. For example, 
as we shall discuss in Sec. 5, the stable primary states of a laser model, for a certain range 
of values of its parameters, are period functions of time [HL1, AS]. 

3. Entropy, Algorithmic Complexity and Disorder. 

Entropy and Disorder. The entropy, S(p), of a state p is given by Von Neumann's 
formula [VN], which, in units for which Boltzmann's constant is log 2 (e), takes the following 
form. 

S(p) = -Tr(plog 2 (p)). (3.1) 
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In order to expose the probabilistic character of this formula, we note that p is a convex 
combination of mutually orthogonal one dimensional projectors, P^, of its eigenvectors, 
i.e. 

Thus w = {wk} is a probability measure on the pure states {Pk} and Eq. (3.1) is equivalent 
to the formula 

S(p) = -^2 k W k log 2 (w k ), (3.3) 

which is just Shannon's formula [SW] for the entropy of the probability measure w. Indeed, 
in the case where E is a classical system and p is a probability measure on a discrete space 
K, the entropy S(p) is given by Eq. (3.3), with Wk the probability attached to the pure 
state represented by the point k of K. 

The formula (3.3), and hence also (3.1), has the natural interpretation [SW, Kh, Sz] 
that —S(p) represents the information carried by the state p; or, equivalently, that the 
value of S(p) is a measure of the degree of disorder of that state. To be precise, it is a 
strictly probabilistic measure of that disorder since , by Eq. (3.3), S(p) depends exclusively 
on the probability measure w and not at all on the structures of the pure states Pj . 

Turning now to the case where E is an infinite system, as formulated in Section 2, 
we denote by p\ the restriction of a state p to the local algebra Aa- The entropy density 
induced by p in the region A is therefore S(pa)/\A\i where the numerator is the Von 
Neumann entropy of p\ and the denominator is the volume of A. It then follows from the 
strong subadditivity of entropy [LR] that, for any translationally invariant state p, this 
local entropy density converges to a limit s(p) as A increases to X over a set of suitably 
regular regions, i.e. 

lim ATX ^^ = S (p) V peS x . (3.4) 

Evidently, it follows from the discussion following Eq. (3.3) that s(p) is a strictly proba- 
bilistic measure of the disorder of the state p. 

Algorithmic Complexity and Disorder. A complementary, intrinsic character- 
isation of disorder, as applied to pure states, has been provided by Kolmogorov [Ko] in 
the classical regime, and quantum versions of this have subsequently been proposed by 
Berthiaume et al [B VL] , Gacs [Ga] and Vitanyi [Vi] . This is based on the concept of algo- 
rithmic complexity, which was introduced by Kolmogorov [Ko] in the context of classical 
communication theory, in the following form. For any string, /cat, of N symbols, drawn 
from the binary set {0, 1}, the algorithmic complexity C(/cjv) is defined to be the length 
of the shortest programme required for the precise specification of that string by a uni- 
versal Turing machine. Such a string of 0's and l's corresponds to a pure state of a one- 
dimensional classical lattice gas, E^r, whose phase space is K N := {0, l}' 1 '^. 

In order to pursue the properties of this algorithmic complexity, we treat Eat as a 
subsystem of the infinitely extended classical lattice gas, E, whose phase space is K := 
{0, 1} Z . Thus the elements of K are the maps k : x^k x of Z into {0, 1} and the spatially 
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ergodic states of E are defined as in Section 2. The following theorem, due to Brudno [Br] 
provides a remarkable relationship between algorithmic complexity and entropy in terms 
of these definitions. 

Theorem 3.1 [Br]. Let p be a spatially ergodic measure on K and, for k<EK , let 
be the restriction of k to [1, AT]. Then, for p-almost all kEK , 

limN^N^CNik) = s(p). (3.5) 



Comment. This thereom signifies that, for p-almost all k in K, the algorithmic 
complexity density C(k]si)/N induced by k on the segment [l,iV] of Z converges to the 
entropy density s(p) as N-^oo. Hence, as algorithmic complexity is an intrinsic measure 
of disorder, the theorem vindicates the standard representation of the disorder of a macro- 
scopic system by its entropy, at least in the case of one-dimensional lattice gases. From the 
physical standpoint, this does not conflict with the fact that the entropy of a pure state is 
zero, for the following reason. The specification of a pure state of a system with N degrees 
of freedom would require the evaluation of N variables. In the case of a macroscopic sys- 
tem, for which N is extremely large (e.g. even in the case of a one-dimensional one, it is 
typically of the order of 10 8 ), such a specification is out of the question. Indeed, for such 
a system, the only accessible information about its state is limited to the determination of 
a 'few'macroscopic variables. The state inferred therefrom is then a highly mixed one. 

Quantum Systems. The above picture of entropy, complexity and disorder has been 
extended to quantum systems in the following way. Firstly, Berthiaume et al [BVL] have 
formulated the algorithmic complexity of a pure state, p, of a string of N qubits as a natural 
quantum analogue of Kolmogorov's classical picture (namely the length of the shortest 
programme required for the determination of p by a universal quantum Turing machine) 
and Benatti et al [BKMSS] have established a quantum version of Brudno's theorem for this 
complexity. These works were formulated within the framework of quantum informatics, 
but they may easily be translated into statistical mechanical terms by noting that a string 
of qubits corresponds to a system of Pauli spins on a one-dimensional lattice. The infinitely 
extended version, E, of such a system is thus a particular case of the quantum lattice model 
of Sec. 2 for which d = 1 and the single particle Hilbert space Ho is two-dimensional. For 
any natural number N, we denote by E^v the subsystem of E comprising the spins at the 
sites 1,2, ,N and by An its algebra of observables. For any state p of E, we denote by pn 
the state of Ejv given by its restriction to An and, following [BKMSS], we term a sequence 
of projectors {PnEAn} p- typical if limjv->oo/o(-FW) = 1. Further, if pn (eAn) is a one- 
dimensional projector we denote by C(j>n) the algorithmic complexity of this state of E^v 
that it represents, as formulated by [BVL]. The quantum version of Brudno's theorem then 
takes the following form [BKMSS]. 

Theorem 3.2. Let p be a spatially ergodic state of the infinite chain, E ; of Pauli 
spins. Then there exists a sequence of p-typical projectors {PnEAn} such that, for any 
e > 0, every one- dimensional projector pn < Pn satisfies the following inequality for N 
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sufficiently large. 



N- 1 C(p N )e(s(p)-e,s(p) + e). 



(3.6) 



Comment. This theorem signifies that the comment following Theorem 3.1 carries 
through to its natural quantum analogue. Further, as pointed out in [BKMSS], the above 
Theorem 3.2 is extendible to any chain of atoms, the observables of each of which are 
irreducibly represented in a finite dimensional Hilbert space Ho- Moreover, we can extend 
that theorem to any <i-dimensional lattice system, E, as formulated in Section 2, in the 
following way. We replace the string [l,iV] of the one- dimensional lattice by the block 
Atv := [1, N] d of Z d in the definitions of E^, Hn and An, so that E^ is now the system 
of atoms occupying the block An- The algorithmic complexity, C(pn) of a pure state pn 
of this system is then the length of the shortest programme required to specify this state 
by a universal quantum Turing machine and its complexity density is C(j>n) /N d . Further, 
if p is a state of the infinite system E, a sequence of projectors {Pn&An} is again termed 
p-typical if limjv_ i . 00 p(Pjv) = l.With these definitions, the treatment of [BKMSS] can be 
carried through to yield the following d-dimensional generalisation of Theorem 3.2. 

Theorem 3.3. Let p be a spatially ergodic state of the infinite system, E, of atoms on 
the lattice Z d . Then there exists a p-typical sequence of projectors {PnEAn} such that, for 
any e > 0, the algorithmic complexity, C(j>n), of any one- dimensional projector pn < Pn 
satisfies the following inequality for N sufficiently large. 



Comment. This theorem signifies that the comments following Theorems 3.1 and 3.2 
extend to quantum systems on lattices of arbitrary finite dimensionality and thus vindicates 
the standard picture wherein the disorder of a state of a quantum lattice system is given 
by its Von Neumann entropy. 



Symmetry Groups and G-fields. A symmetry group of the model E = (A, 5, 7, a) 

is a group, G, that has a faithful representation, 9, in Aut(^4). It is a dynamical symmetry 
group if 9(G) commutes with the time-translations at- For any symmetry group G, an 
n-component quantum field £ = (£1, . ., £ n ), affiliated to A, is termed a G-field if the action 
of 9(G) on £ takes one of the following forms. 

(a) In the case where this symmetry group is spatial, 



N- d C(p N )e(s(p)-e,s(p) + e). 



(3.7) 



4. Symmetry and Order. 



[6(g)t](x) = t(T g x) VgeG, 



(4.1a) 



where T g is a transformation of X. 
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(b) In the case where G is an internal symmetry group, 

[6(s)t]j(*) = T, n k=1 V 9^k(x) V geG, (4.16) 

where V g = [V g jk] is a unitary transformation of R n or C n , according to whether the field 
£ is real or complex.. 

Thus, in either case, the action of 6(g) on £ is that of a linear transformation <p g of this 
field alone, without involvement of other observables of the system, i.e. 

[0(9)t](*)=[M(x) •■=&*)■ (4-1) 

Symmetry Breakdown. An stationary state p of S that is not invariant with respect 
to a dynamical symmetry group G is said to spontaneously break that symmetry. In this 
case, the states {p g := po0(g)\g£G} are also G-symmetry breaking stationary states of the 
system. We term this set of states the G-orbit of p and denote it by Oq(p)- 

In the case where E is conservative and p&£(Sp), the set of extremal KMS states at 
inverse temperature f3, it follows from the KMS condition (2.1) and the definition of Oq{p) 
that £(S(3)<zOg{p)- Further, defining 

l(x) := (p;£(:r)) and := (p ff ; £(*)), (4.2) 

it follows from Eq. (4.1) that 

where the action of <p g on the classical field £ is the same as that on the quantum field 
£. We assume that this action is non-trivial* unless g is the identity element of G and 
consequently that £ g ^£ g > if g^g' ■ Hence, by Eqs. (4.1)-(4.3), the G-field £ serves to 
separate the states of the orbit Oq(p)- We note here that, since these states are primary, 
it follows from Eqs. (2.4) and (4.1)-(4.3) that 

lim a ^ 00 [(p; %(x)4 g (x + a)) - T g (x).£ g (x + a)] = 0, (4.4) 

where the dot denotes the scalar product in R n or C n . 

Symmetry Breakdown and Order. At a qualitative level, we conceive a state of 
a complex system to be ordered if its components cooperate in such a way as to produce a 
macroscopic field or signal. Thus, a prototype example of ordering in this sense is that of 
the alignment of the spins (mini-magnets!) of a ferromagnetic material so as to produce a 
resultant polarisation. Hence, assuming that the dynamics of the material is rotationally 

* In fact, this condition is effectively fulfilled if the subgroup, H, of G under which the 
field £ is invariant is a normal one. For in that case we may replace the group G of the 
present treatment by the factor group G/H, which then satisfies our demands. 
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invariant, this magnetic ordering amounts to a breakdown of its rotational symmetry, as 
manifested by the direction of the polarisation field. 

The generalisation of this concept of order through symmetry breakdown is straight- 
forward. Thus, an equilibrium state p of £ is said to be ordered if it is not invariant with 
respect to a dynamical symmetry group G, and its ordering is then represented by the 
^-dependence of the classical field £ . Among the numerous proven examples of such order 
are the following. 

Example 1. The two-dimensional Ising model, whose pure equilibrium phases at any tem- 
perature below its critical point are polarised [MM] . For this model, G is the binary group 
(e,r), where e is the identity element and 9(r) is the spin reversal automorphism. The 
G-field is the Ising spin s x at the site x. 

Example 2. The antiferromagnetic phase of the Heisenberg model [DLS]. For this, G is 
the three dimensional rotation group and 9(G) represents its action on the Pauli spins 
constituting the model. The G field is the Pauli spin vector a(x) at the site x and so 
£ g (x) is a spatially periodic vector filed whose direction is determined by g. 

Example 3. A crystalline phase of the generic model of a continuous system £ that occupies 
the space X = H d . For this system, G is the factor group X/Y, where X is the additive 
space translation group and Y is the normal subgroup of X corresponding to the crystal 
structure [EKV]. The G-field £(x) may be chosen to be the particle density ip*(x)ip(x). 

Example 4- This is the order corresponding to a generalised version of Bose-Einstein (BE) 
condensation, which was first proposed by O. Penrose and Onsager [PO] as a characterisa- 
tion of the superfluidityof Hell and subsequently extended by Yang [Ya] to superconduc- 
tors. In fact, this order corresponds to the breakdown of gauge symmetry*, the relevant dy- 
namical symmetry group, G, being that of the gauge transformations tp(x)— ^ip(x)(exp(ic), 
where c runs through the reals. In the case of bosons, the G-field is just ip(x): in the 
case of fermions it is the pair field ip^(x)ipi(x), the two factors being the components of 
ip(x) with spin parallel and antiparallel to some fixed axis. 

Example 5. This is the pumped phonon model, which was proposed by H. Froehlich [Fr] in 
a biological context and put onto a rigorous footing by Duffield [Du] . It is an open dissative 
system, consisting of N phonon modes that are coupled to energy pumps and sinks and 
that exchange quanta with one another in conformity with the principle of detailed balance 
at the prevailing environmental temperature. Remarkably, the model is driven by these 
forces into a state wherein a macroscopic number, of order N, of its quanta condense into 
the mode of lowest frequency when the pumping strength exceeds a critical value. Thus, 
in the limit A^oo, it exhibits a BE condensation into a nonequilibrium steady state. This 
amounts to a gauge symmetry breakdown, far from thermal equilibrium. The source of this 
phenomenon stems is the competition between the pumping and the discharge of quanta 
into the sinks, which fixes the total number of quanta and thereby renders the state of the 

* See [Sel, Ch.9] for a detailed discussion of this symmetry breakdown and its connec- 
tions with both superfluidity and the so-called 'off-diagonal long range order' of O. Penrose 
and Onsager [OP] and Yang [Ya]. 
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model similar to an equilibrium state of an ideal Bose gas with a fixed number of particles. 
By contrast, a phonon system at equilibrium with a thermostat at fixed temperature has 
a Planck distribution of its quanta and thus does not experience a BE condensation. 

5. Coherence. 

Coherence is an extreme version of order. Specifically, a state p of the system E is 
said to be coherent with repect to a G-field, if it satisfies Glauber's [Gl] condition 

that 

<p;£ # (*i). .i*{x n )) = U] =1 (p : £#(*,■)>, (5-1) 

where each is either £ or £* and (p; is a non-trivial function of x, which is simple 
in that it involves just a 'few' parameters: the simplicity condition serves to exclude cases 
where the classical field represented by this function varies chaotically with x. 

The natural dynamical version of this coherence condition is that obtained by the 
replacement, in Eq. (5.1), of by its evolute (= ct t ^{x)). Thus, the resultant 

dynamical coherence codition is that 

(p;&f (*0- •££(*»)> = n?=i<p = (5-2) 

where (p; £t(x)) is a non-trivial, simple function of both x and t. 

It follows from these specifications that coherence corresponds to the behaviour of 
the quantum field or £t(x), in the state p, as a classical, dispersion-free field := 

(p;£,{x)) or f t (x) := (p;£,t{x)) 

The following examples of coherence have been established in different variants of the 
Dicke model [Di] of matter interacting with a single radiative mode. 

Example 1. Hepp and Lieb [HL2] have proved that the Dicke model has a low tem- 
perature equilibrium phase characterised by super-radiance, i.e. by coherence of the pure 
equilibrium phase with respect to the radiation field, with breakdown of the gauge sym- 
metry of that field. They also treated [HL1] the open version of the Dicke model in which 
each atom of the matter was coupled to an energy pump and a sink and the radiation 
mode was coupled to a sink. This model undegoes a phase transition, far from thermal 
equilibrium, at a critical value, p c , of the pumping strength, p. Specifically, for p < p c , its 
stable state is stationary and its radiation incoherent; while for p > p c , the stable state is 
simply periodic in time and its radiative mode coherent, again with breakdown of gauge 
symmetry. Thus the transition is from normal light to laser light. 

Example 2. Alii and Sewell [AS] extended the open version of the Dicke model to 
one with many modes, each with its own sink, and obtained the following picture of its 
phase structure in terms two critical values, p\ and pi (> pi) of the pumping strength. 
For p < pi, there is a unique stable stationary state of the model, and the radiation is 
normal, i.e. incoherent. For pi < p < p 2 , the stable state varies periodically with time and 
the radiation is coherent and monochromatic, again with breakdown of gauge symmetry. 
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For p > P2, the radiation is chaotic according to either the mechanism of Ruelle- Takens 
[RT] (strange attractor) or that of Landau [LL] (multimode turbulence), depending on the 
parameters of the model. 

Example 3. Pule, Verbeure and Zagrebnov [PVZ] constructed a model of a system 
of interacting two level bosonic atoms that is coupled to a radiative mode. This coupling 
was shown to leads to a rich equilibrium phase structure. In particular, for sufficiently 
large values of the chemical potential, the radiation is coherent and the matter exhibits a 
two-fold BE condensation, one for each of its atomic levels. This supports the experimental 
observation that the action of a laser field on a bosonic condensate of atoms with internal 
structure leads to enhancements of both the laser field and the BE condensation [KI] . 

6. Concluding Remarks 

We may summarise the contents of this article as follows. 

Disorder corresponds to randomness, of which entropy and algorithmic complexity 
provide probabilistic and intrinsic measures, respectively. Remarkably, Brudno's theorem 
and its quantum analogues have established the essential equivalence of these measures in 
the thermodynamic limit (cf. Theorems 3.1-3.3). 

Order, on the other hand, constitutes organisation manifested by a macroscopic field or 
signal, and this can prevail amidst high disorder. Thus, in the present physical context, the 
concepts of order and disorder are certainly not antitheses of one another. The particular 
types of order and coherence described here stem from symmetry breakdown and BE 
condensation, and certainly do not cover all the kinds of organisation that arise in complex 
systems. For example, it is still true that, as pointed out by Schroedinger many years ago 
[Sc], the presently available pictures of order do not cover biological organisation.. Thus, a 
major challenge of statistical physics is to characterise other kinds of ordering that exist 
in nature. 

Acknowledgement. It is a pleasure to thank Fabio Benatti for an enlightening 
correspondence on the subject of algorithmic complexity. 
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